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^ : ABSTRACT 

We have discovered two simultaneous kHz quasi-periodic oscillations 
(QPOs) in the Z source GX 340+0 with the Rossi X-ray Timing Ex- 
! plorer. The X-ray hardness-intensity and color-color diagram each show 
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a full Z-track, with an extra limb branching off the flaring branch of 
the Z. Both peaks moved to higher frequencies when the mass accretion 
rate increased. The two peaks moved from 247 ± 6 and 567 ± 39 Hz at 
the left end of the horizontal branch to 625 ±18 and 820 ± 19 Hz at its 
right end. The higher frequency peak's rms amplitude (5-60 keV) and 
FWHM decreased from ~ 5% and 383 ± 135 Hz to ~ 2%, and 145 ± 62 
Hz, respectively. The rms amplitude and FWHM of the lower peak were 
consistent with being constant near 2.5 % and 100 Hz. The kHz QPO 
separation was consistent with being constant at 325 ± 10 Hz. 

Simultaneous with the kHz QPOs we detected the horizontal branch 
oscillations (HBO) and its second harmonic, at frequencies between 20 
and 50 Hz, and 38 and 69 Hz, respectively. The normal branch oscilla- 
tions were only detected on the upper and middle normal branch, and 
became undetectable on the lower normal branch. The HBO frequencies 
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do not fall within the range predicted for Lense-Thirring (LT) precession, 
unless either the ratio of the neutron star moment of inertia to neutron 
star mass is at least 4, 10 45 e/ cm 2 /M Q , the frequencies of the HBO are in 
fact the sub-harmonic oscillations, or the observed kHz peak difference 
is half the spin frequency and not the spin frequency. 

During a 1.2 day gap between two observations, the Z-track in the 
hardness-intensity diagram moved to higher count rates by about 3.5%. 
Comparing data before and after this shift, we find that the HBO prop- 
erties are determined by position on the Z-track and not directly by 
count rate or X-ray colors. 

Subject headings: accretion, accretion disks - stars: individual (GX 340+0) 
- stars: neutron - X-rays: stars 
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1. Introduction 

GX 340+0 is a bright low-mass X-ray bi- 
nary and a Z source flHasinger fc van der 
Klis 1989] ). The Z-track traced out in the X 



ray color-color (CD) and hardness-intensity 
diagram (HID) is usually built up of three 
components, the horizontal branch (HB), the 
normal branch (NB), and the flaring branch 
(FB). The mass accretion rate (M) is thought 
to increase from a few tens of the Edding- 
ton mass accretion rate when the source is on 
the HB, to near-Eddington on the NB, and 
super-Eddington on the FB. On the HB and 
the upper part of the NB of GX 340+0, quasi- 
periodic oscillations (QPOs) with frequencies 
of 32-50 Hz, the horizontal branch oscilla- 
tions (HBOs), were detected by Penninx et 
al. (1991). On the central part of the NB, 
QPOs with a typical frequency of 5.6 Hz, the 
normal branch oscillations (NBOs), were de- 
tected by van Paradijs et al. (1988). 

In four Z sources, ScoX-1 ( [van der Klis et 
al. 1996at |van der Klis et al. 1997|) , GX5-1 



van der Klis et al. 1996b| ), GX17+2 (|WTp 



aands et al. 1997b| ), and CygX-2 ([Wijnands 
|et al. 1998a|) two simultaneous kilo Hertz 
(kHz) QPOs have been observed, with fre- 
quencies reflecting flux changes on timescales 
close to the dynamical timescales near a neu- 
tron star. The higher frequency peak (hence- 
forth called upper peak) has been proposed to 
reflect the Keplerian frequency of blobs of ma- 
terial orbiting around the neutron star, at a 
preferred radius in the disk, such as the sonic 
radius (Miller, Lamb & Psaltis 1998). The 
lower frequency peak (henceforth called lower 
peak) is the proposed beat frequency between 
the frequency of the upper peak and the spin 
frequency of the neutron star. The kHz QPO 
frequency is observed to increase with increas- 
ing M. Similar kHz QPOs are also found 



in atoll sources. The frequencies of the kHz 
QPOs are between 300 and 1200 Hz. (See van 
der Klis (1997) for a recent overview). 

In this Letter we report the discovery of 
two simultaneous kHz QPOs in the Z source 
GX 340+0. A preliminary announcement of 
these results was already made by Wijnands 
& van der Klis (1998b). 

2. Observations and Analysis 

We observed GX 340+0 with the propor- 
tional counter array onboard the Rossi X-ray 
Timing Explorer on 1997 Apr. 17, Sept. 21, 23, 25, 
and Nov. 1, 2, 3, and 4. The total amount of 
good data obtained was 138 ksec. During 5% 
of this time only three or four of the five de- 
tectors were active. For our CD and HID we 
used only the data from the three detectors 
that were always active. In our power spec- 
tral analysis we used all available data. The 
mean five-detector countrate was 6443 cts s _1 . 
The data were obtained using 16s time reso- 
lution in 129 photon energy bands (effective 
energy range 2-60 keV), and simultaneously 
with 122 fj,s time resolution in four bands (2- 
5.0 keV, 5.0-6.4 keV, 6.4-8.6 keV, and 8.6-60 
keV) on Apr. 17 or in two bands for the other 
dates (2-5.0 keV and 5.0-60 keV). 

We calculated power density spectra us- 
ing the 122 fxs data divided into 16s seg- 
ments. To determine the properties of the 
kHz QPOs we fitted the 95-4096 Hz power 
spectra with a function that consisted of one 
or two Lorentzians (the kHz QPOs) and a 
constant plus a broad sinusoid (the dead-time 
modified Poisson noise) ( |Zhang et al. 1995 ). 



The Very Large Event window (|van der Klis 
et al. 199TD was set to 55 /is, so that its effect 
on the Poisson noise is small and could be in- 
corporated in the sinusoid. To determine the 
properties of the HBO, we fitted the 0.125- 
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126.5 Hz power spectra using a function which 
consisted of the sum of a constant, one or two 
Lorentzian peaks, and two exponentially cut- 
off power-law components. The errors were 
determined using A% 2 = 1.0 (la single pa- 
rameter) and upper limits using A% 2 = 2.71, 
corresponding to a 95% confidence level. Up- 
per limits on the kHz QPOs were determined 
using a fixed FWHM of 150 Hz. When only 
the lower peak was detected, the upper limit 
on the upper peak was determined by setting 
the frequency to the value of the frequency of 
the lower peak plus the mean frequency dif- 
ference between the two peaks. Upper limits 
on the HBO were determined by keeping the 
FWHM fixed at 10 Hz. 

In the HID and CD, the hard color is de- 
fined as the logarithm of the 9.7-16.0 keV/ 6.4- 
9.7 keV countrate ratio. The soft color is de- 
fined as the logarithm of the 3.5-6.4 keV/2- 
3.5 keV countrate ratio. The intensity is de- 
fined as the logarithm of the countrate in 
the 2-16.0 keV band. The HID and CD are 
background-corrected but no dead-time cor- 
rection was carried out. The dead-time cor- 
rection factor was ~ 2%. In order to measure 
the position along the Z-track, we use the S z 
parametrization (e.g. Wijnands et al. 1997a ) 
applied to the HID. In this parametrization 
the HB/NB vertex is fixed at S z = 1.0 and 
the NB/FB vertex at S z = 2.0. We selected 
the power spectra according to S z , and deter- 
mined the average S z and standard deviation 
for each selection. 

The track of the Nov. 4 data is shifted to- 
wards higher intensity and softer colors with 
respect to the Z-track in the other data. We 
therefore determined the HBO properties for 
the Nov. 4 data (~ 18 ksec) separately. To 
determine the kHz QPO properties we com- 
bined the Nov. 4 data with the rest in order to 
get a higher signal to noise ratio. As the kHz 



QPOs were only significantly detected in the 
5.0-60 keV band and the significance did not 
increase combining the two bands, we used 
only the 5.0-60 keV band in our analysis. 

3. Results 

We observed a full Z-track in the HID and 
CD (Fig. |l|). The Z-shape shows a fourth 
branch beyond the FB. This extra, trailing 



branch was seen before in GX 340+0 (|Fem 
ninx et al. 1991| ) and GX5-1 ( |Kuulkers et al 



19941 |Wijnands fc van der Klis 1998b|) . De- 



tailed examination of the HB/NB vertex in 
the HID revealed that the 1997 Nov. 4 data 
were shifted to higher count rates (~ 3.5%) 
and softer colors (~ 0.5%) than the other 
data. Due to this shift, which occurred dur- 
ing a 1.2 day gap between two observations, 
the HB/NB vertex in the HID is broadened. 
Due to the size of the dots in Fig. |Tj this is not 
visible. This shift was also seen in EXOSAT 
data of GX 340+0 ( |Kuulkers & van der Klis 



19961) . 



We detected kHz QPOs in the 5.0-60 keV 
band (Fig. 0). The lower kHz peak was de- 
tected at S z values up to 1.04, the upper peak 
up to S z = 0.95 (Fig. [3]). The lower peak had 
a frequency of 247 ± 10 Hz at S z = 0.47 and 
increased to 625 ± 18 Hz at S z = 1.04, its rms 
amplitude remained about constant at values 
near 2%, and its FWHM varied between 30 
and 160 Hz with no clear correlation with S z . 
The frequency of the upper peak was 567 + 39 
Hz at S z = 0.47 and increased to 820 ± 19 
Hz at S z = 0.95, while its rms amplitude de- 
creased from 5.1 + 0.9% to 2.0 + 0.4%, and its 
FWHM from 383+135 Hz to 145+62 Hz. The 
FWHM ratio did not significantly depend on 
S z . When both peaks were detected, the peak 
separation was consistent with being constant 
at a value of 325 ± 10 Hz. We did not de- 
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tect significant kHz QPOs in the 2-5.0 keV 
bands, with upper limits of 2.8% and 2.3% at 
S z = 0.58 and 2.7% and 2.8% at S z = 0.87 
for the lower and upper peaks, respectively, 
while in the 5.0-60 keV band the peaks were 
found with a rms amplitude of 1.5% and 4.4% 
at S z = 0.58 and 1.6% and 2.8% at S z = 0.87, 
respectively. 

Simultaneously with the kHz QPO, the 
HBO and its second harmonic were detected. 
When we used the same S z selections as in 
the fitting of the kHz QPOs the HBO some- 
times showed double peaks, especially at low 
S z . This was caused by the steep depen- 
dence of the HBO frequency on S z and its 
narrow peak width. Therefore up to S z = 0.7, 
we used narrower S z selections in fitting the 
HBO to prevent the HBO peak from mov- 
ing too much within the selection. The HBO 
frequency increased from 19.43 ± 0.03 Hz at 
S z = 0.43 to 49.92 ± 0.21 Hz at S z = 1.05. 
Then it remained constant up to S z = 1.25 
(Fig. |3] b). The rms amplitude in the 5.0- 
60 keV energy band of the HBO decreased 
smoothly from 9.1 ± 0.1% at the left side of 
the HB (S z = 0.43) to 1.8 ± 0.2% on the NB 
(S z = 1.25). For higher S z selections, no HBO 
was detected, with upper limits of ~ 2% rms. 
The FWHM increased from 2.9 ± 0.1 Hz to 
15.3 ±4.2 Hz. We did not find a clear depen- 
dence on S z of the rms amplitude and FWHM 
of the second harmonic. Due to the presence 
of a broad noise component around the fre- 
quency of the second harmonic it was diffi- 
cult to determine the properties of either one 
of the two components. 

We found that the HBO frequency vs. S z 
relation of the Nov. 4 data was offset with re- 
spect to that of the other data (Fig. |] left), 
due to the shift in the vertex of the Nov. 4 
data, which along the HB, corresponds to a 
shift in S z of ~ 0.05, equal to the offset in 



the S z vs. vubo plot (Fig- ||) left). If we 
correct for the change in S z by measuring S z 
on Nov. 4 from the HB/NB vertex appropri- 
ate for that date, the HBO frequencies, rms 
amplitude, and FWHM of the Nov. 4 data 
and the HBO frequencies, rms amplitude, and 
FWHM of the other data as a function of S z 
are the same (Fig. |], right). This strengthens 
previous conclusions flKuulkers et al. 1994p 
that S z , not count rate, is the better measure 
for M. 

We also detected the NBOs, between S z ~ 
1.1 and S z ~ 1.8. These will be more exten- 
sively discussed elsewhere. 

4. Discussion 



We have discovered kHz QPOs in the Z 
source GX 340+0. The frequencies of both 
the lower and upper peak increased as the 
source moved along the HB to the HB/NB 
vertex. The peak separation was consistent 
with being constant. This is similar to what 
was found in GX 17+2 (|Wijnands et al. 1997b|) , 



CygX-2 ([Wi.inands et al. 1998a| ), and GX5-1 



( |van der Klis et al. 1996b| ; |Wijnands et al 



[T998"cD . Only in ScoX-1 (|van der Klis et al 



1997| ) and perhaps in the atoll source 4U 1608- 



52 ( Mendez et al. 1998|) the peak separation 
was found to decrease with increasing M . If 
this decrease in Sco X-l is in some way related 
to the source approching the Eddington criti- 
cal luminosity ( White fc Zhang 1997| ; [Miller et 
al. 1998|) , this could mean that if we were able 



to detect the kHz QPOs in GX 340+0 and the 
other Z-sources (further) up the NB we might 
also see the peak separation decrease. 

The lower peak reaches the so far low- 
est "kHz" QPO frequency found in any low 
mass X-ray binary, 247 Hz. The maximum 
frequency reached by the upper peak, 820 
Hz, is also rather low. Using the sonic-point 
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model flMiller et al. 1998|) to explain the 
upper peak kHz QPO frequency, we found 



lo\ter limits on the mass of the neutron star 



of M neutr > O.6M for a neutron star radius 
of 4M, and M neutr > 1.2M Q for a neutron 
star radius of 7M. (The implicit assumptions 
made in calculating these lower limits are: 
Schwarzschild spacetime metric, non-rotating 
radiating layer, and photons remove angular 
momentum from the disk material only in one 



sc attering event.) We note that the frequem 



cies of the kHz QPOs seem to confirm the idea 
(e.g. IKuulkers et al. 1994|) that GX 340+0 is 
more similar in appearance to GX5-1 ( |Wij- 



na hds et al. 1998c]) , than to ScoX-1 (|van der 



Kl |s et al. 1991 ), and GX 17+2 flWijnands et| 
al. 1 1997b]) (frequencies higher than 1000 Hz). 



Any straightforward beat-frequency model 



predicts the upper peak to be narrower than, 
or as narrow as the lower peak. We find 



in |GX 340+0 that the upper peak is usually" 



broader than the lower peak. In fact, the up- 
per peaks found in GX 340+0 were very broad 
with up to 0.67 at S z = 0.43. Therefore a 
straightforward beat-frequency interpretation 
of the two peaks in GX 340+0 is difficult. In 
GX 340+0, GX5-1 QWijnands et al. l998c| ), 
and GX17+2 QWijnands et al. 1997b| ) the 
upper peak FWHM decreases as a function of 
S g . To explain this behaviour one could in- 
volve scattering in a rapidly variable medium 
( |Mendez et al. 1998| ). The effect of this scat- 
tering should become less as M increases. 

Although the behaviour of the kHz QPOs 
as a function of inferred M in GX 340+0 
is similar to that found in other Z sources, 
the HBO behaviour is rather different. The 
FWHM of the HBO increased as the rms 
amplitude decreased, which is unlike what 
is seen in other Z sources. In CygX-2 Wi- 
jnands et al. (1997c) found the FWHM to 
be well-correlated to the rms. In GX 17+2 



the FWHM remained about constant while 
the rms amplitude decreased (Homan et al 



1998| ). In previous observations of GX 340+0 



( Penninx et al. 199 1| , Table 2) the behaviour 
is consistent with our findings. In obser- 
vations of Kuulkers & van der Klis (1996) 
the FWHM does not show a clear relation 
as a funtion of the rms amplitude. Ex- 
actly the same behaviour (as we reported for 
GX 340+0) was found in GX5-1 ((Lewin et al 



19921 IWijnands et al. 1998c|) . This suggests 



that an extra parameter in addition to M de- 
termines the timing properties. Another indi- 
cation of this is that there are differences be- 
tween the Z sources in the position on the Z- 
track where the NBO are found. In GX 17+2 
( Penninx et al. 1990Q , and ScoX-1 ( |van dei 
Klis et al. 1987| ) the NBO were detected up 
the FB, e.g. S z > 2. In CygX-2 they were 



detected up to the NB/FB vertex ( |Wijnandi 



et al. 1997a| ), while in GX 340+0 and GX5-1 
the NBO were detected up to S z = 1.8. 

Recently it has been proposed that the 
HBO frequency may reflect the general rel- 
ativistic Lense-Thirring (LT) precession fre- 
quency at the inner edge of a warped accre- 
tion disk ( Stella fc Vietri 1998| ). If the up- 
per peak is interpreted as a Keplerian fre- 
quency in the inner part of the disk, there is a 
quadratic relation between the LT precession 
frequency and the upper peak frequency. Ap- 
plying this theory to our observations and as- 
suming the peak difference to be the neutron 
star spin frequency, the LT predictions for the 
HBO frequency are too low by at least a fac- 
tor of 2. The same discrepancy by a factor of 
2 has been reported for the Z source GX 17+2 
( (Stella & Vietri 1998| ) . There are several ways 
out, but each appears to create serious diffi- 
culties. I/M (with I the moment of inertia 
of the neutron star and M its mass) may be 
a factor of 2 higher than presently thought, 
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however, no sensible equation of state appears 
to be able to produce such values of I/M. 
The kHz QPO peak difference may not be 
equal to the neutron star spin frequency but 
to half the spin frequency, however, at least in 
atoll sources this seems to contradict the re- 
sults obtained from burst QPO ( |Strohmayer 
et |al. 1996|) . The HBO may not be the fun- 



damental frequency but the harmonic of the 
predicted LT frequency, however, in none of 
the Z sources there is any evidence of excess 
power at frequencies half the HBO frequen- 
cies. 
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Fig. 1. — Hardness-intensity (left) and color- 
color diagram (right) of GX 340+0. The hard 
color is defined as the logarithm of the 9.7- 
16.0 keV / 6.4-9.7 keV count rate ratio and 
the soft color as the logarithm of the 3.5-6.4 
keV/ 2-3.5 keV countrate ratio. The three- 
detector countrate is measured in the 2-16.0 
keV band. The data was background sub- 
tracted but no dead-time correction was ap- 
plied. The typical error bars are visible in the 
HID and CD in points in the extra trailing 
branch. 

Fig. 2. — Typical power density spectra of 
GX 340+0; the power is in units of fractional 
amplitude squared per Hertz. The kHz QPOs 
are at frequencies of 452 and 753 Hz, respec- 
tively, at S z = 0.87. The rise in the power 
density of the left plot towards higher fre- 
quencies is caused by instrumental dead-time 
effects on the Poisson noise. 

Fig. 3.— (A) Frequencies of the kHz QPOs, 

(B) frequencies of the HBO and its harmonic, 

(C) rms amplitude of the upper kHz QPO, 

(D) rms amplitude of the HBO, (E) rms 
amplitude of the lower kHz QPO, and (F) 
FWHM of the HBO, as a function of S z . The 
error bars on S z are standard deviations. 

Fig. 4. — Left: The HBO frequencies of the 
November 4 data (circles) and the other data 
(X-es). Right: The frequencies of the HBO 
for the November 4 data corrected for the shift 
(circles) and the other data (X-es). The error 
bars in the frequency are smaller than the size 
of the symbols, (see also Section 3 paragraph 
4). 
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